Periodically modulated quantum 
nonlinear oscillators 

M. I. Dykman 



7.1 Introduction 



Vibrational systems have been attracting much attention in physics. Such systems are 
always nonlinear, at least to some extent. For weak damping, even small nonlinearity 
can become important. For example, classically, the nonlincarity-induccd dependence 
of the vibration frequency on amplitude can lead to bistability of forced resonant vi- 
brations ( [Landau and Lifshitz, 2004] ), see Fig. 17.11 Quantum mechanically, the nonlin- 
earity makes the frequencies of transitions between adjacent energy levels different and 
thus enables spectroscopic identification and selective excitation of these transitions. 
The interest in quantum effects in oscillators significantly increased recently in the 
context of nonlinear vibrations in Josephson junction based systems and applications 



of these systems in quantum information (Wallraff et ai, 2004 Siddiqi et ai, 2005 


Lupa§cu et al, 2006 Steffen et ai, 2006 


Metcalfe et ai, 20071 |Schreier et al, 20081 |Watanabe et ai, 2009| 


Mallet et al., 2009; Vijay et ai, 2009; 'Wi 


Ison et al, 2010t|Bishop et ai, 2010 Reed et al, 2010). 



reached also in nanomechanical resonators (O'Connell et ai, 2010 Riviere et ai, 2011) 



This development makes it possible to study quantum effects in individual vibrational 
systems rather than ensembles. 

Besides being interesting on their own, nonlinear oscillators allow addressing some 
fairly general physics problems. One of them is classical and quantum fluctuations 
far from thermal equilibrium and whether they have features that have no analog in 
systems close to equilibrium. A resonantly modulated nonlinear oscillator could be the 
first well-characterized physical system with no detailed balance, which was used to 
study an important class of fluctuation phenomena, the fluctuation-induced switching 
between coexisting stable states, both in the classical and quantum regimes, and to 
reveal some of such feature^ 



^For the theory and experiment on switching of a resonantly modulated oscillator with no de- 
tailed balance see, in pa rticular, [Dykman and Krivoglaz 1979] [Dmitriev and Dyakonov 1986a 
Vogel and Risken 1990] [Dykman ef aL 1998| [Lapidus ef a/. 1999| [Siddiqi et 2005 



Kim et al. 2 005' 'Aldridgc and Cleland 2005[ 



classical and Dykman and Smelyansky 1988 
Marthaler"a nd Dykm an 2006| 
Mallet et al. 2009, " 



Stambaugh and Chan 2 006a Almog ei al. 2007 for the 
[Vogel and Risken 198 8 Kinslcr and Drumm ond 19^1 



IKatz et al. 20071 ISerba n and W ilhelm 20071 

Pcano and Thorwart 2010al IWilson et'al. 2010. for the quantum regime. 



[Vijay et al. 2009 
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An important source of quantum fluctuations in an oscillator is coupling to a 
thermal bath. The coupling leads to oscillator relaxation via emission of excitations 
in the bath (photons, phonons, etc) accompanied by transitions between the oscillator 
energy levels. If the coupling is weak, the transition rates are small compared to the 
energy transferred in a transition, in frequency units. In the classical limit, the energy 
levels are not resolved and the transitions lead to friction. 

At the quantum level, one should take into account that the transitions happen 
at random. The randomness gives rise to a peculiar quantum noise. For a resonantly 
modulated oscillator, it leads to diffusion over the quantum states, which are time- 
dependent because of the modulation. The result of this diffusion is quantum heating 
of the oscillator. It is seen, in particular, in a nonzero width of the distribution over the 
oscillator states even where the temperature of the thermal reservoir is zero. Quantum 
heating is qualitatively different from the familiar Joule heating, which characterizes 
the power absorbed from the modulating field. In contrast to the Joule heating, the 
resulting distribution over the states turns out to be independent of the oscillator 
relaxation rate, for weak damping. 
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Fig. 7.1 Bistability of forced vibrations for additive (a) and parametric (b) modulation at 
frequency top close to the oscillator eigenfrequency ujo and 2ajo, respectively; Aosc is the scaled 
vibration amplitude. In (a), /3 is the scaled squared modulation amplitude, and in (b) fj, is the 
scaled frequency detuning ujf ~ 2ujo. Parameters /i and the scaling factor of the vibration 
amplitude C are defined in Table 17.11 Solid and dashed lines show the stable and unstable 
stationary vibrational states; in (b) there are two states with the same nonzero Aosc and the 
phases that difi'er by tt. The bifurcation parameter values /3bi,2, Msi,2 indicate where the 
number of stable vibrational states changes. The scaled decay rate k in (a) and (b) is 0.25 
and 0.3, respectively. 



A consequence of quantum heating is quantum activation ( Dykman and Smelyans ky, 1988[ 



[Marthaler and Dykman, 2006| [Dykman, 20071 |Katz et ai, 2007| [Peano and T horwa rt, 2010&D 



This is a mechanism of switching between coexisting stable states of forced vibrations 
shown in Fig. 17.11 Similar to quantum heating, the switching is due to the quantum 
noise that accompanies relaxation. It occurs via transitions over an effective barrier 
that separates the vibrational states. The mechanism differs from quantum tunneling, 
and moreover, leads to an exponentially larger switching rate even for low tempera- 
tures. It differs also from thermal activation, which becomes important in the classical 
regime of high temperatures. 

Quantum activation has no analog in systems close to thermal equilibrium. The 
switching rate has a characteristic dependence on the parameters and displays scaling 
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behavior with characteristic exponents. This has made it possible to identify quantum 
activation in the experiment ( [Vijay et al., 2009] ) . The physics of quantum heating and 
quantum activation is explained in Sees. W72\ and [7731 respectively. 

Quantum heating is manifested also in the power spectra of resonantly modulated 
oscillators and the spectra of their response to a weak additional field. Spectroscopy 
has been recognized as a means of studying the dynamics of modulated oscillators and, 



more recently, of using oscillators for quantum measurements ( Dykman and Krivoglaz, 1979| 



Drummond and Walls, 1980 Collett and Walls, 1985 Dykman et al, 1994 Stambaugh and Chan, 20066 



Chan and Stambaugh, 2006 Nation et al, 2008 Wilson et ai, 2010 Vierheilig and Grifoni, 2010 



Boissonneault et ai, 2010) Lafiamme and Clerk, 2011 1. Because of the interplay of 



many interstate transitions at close frequencies, the spectra of modulated quantum os- 
cillators have a characteristic shape. They can have a fine structure, in which case they 
directly provide the quantum-heating induced distribution over the oscillator states 
( [Dykman et ai, 201 1[ ) . The power spectra are also important for understanding the dy- 
namics of resonantly driven oscillators coupled to two-level systems ( [Picot et ai, 20081 
[Serban et al., 2010[ ); interesting spectral manifestations of quantum heating in such 
coupled systems have been recently founcH. The spectra of modulated oscillators are 
discussed in Sec. 17.41 

Nonresonant modulation can also have pronounced effect on the oscillator dy- 
namics. Recently such modulation attracted much attention in optomechanics, where 
intracavity modes are coupled to mechanical vibrations, for example, to the vibrations 
of a mirror in the cavity ( [Kippenberg and Vahala, 2008[ ). Modulation can lead to cool- 
ing and heating of an oscillator, or excite self-sustained vibrations. In contrast to the 
quantum heating discussed above, where quantum fluctuations broaden the distribu- 
tion over the oscillator states in a strong resonant field, here the issue is the change of 
the distribution over the Fock states of the oscillator. An interesting feature of heating 
and cooling is that, in the important case where the effective friction is linear (as for 
standard viscous friction), the distribution over the Fock states remains of the Boltz- 
mann form, but the temperature differs from the bath temperature ( [Dykman, 1978 



Clerk, 2004[ [Wilson-Rae et ai, 2007) [Marquardt et ai, 2007[ ). In Sec. [LS we extend 



the previous analysis (Dykman, 1978) to take into consideration both the nonlinearity 
of the coupling of the oscillator to other degrees of freedom and the nonlinearity of 
the interaction with the modulating field. 



7.2 Resonant modulation: Quantum heating 

For weakly damped oscillators, vibration nonlinearity becomes important once the 
change of the vibration frequency due to the nonlinearity Auj becomes comparable to 
the oscillator decay rate F, which characterizes frequency uncertainty. This happens 
where Aw is still small compared to the oscillator eigenfrequency wq and the nonlinear 
part of the vibration energy is small compared to the harmonic part. Respectively, the 
vibrations remain almost sinusoidal, which significantly simplifies the analysis. At the 
same time, for weak damping an already moderately strong resonant modulation can 



^F. R. Ong et al., in preparation (experiment) and M. Boissonneault et al., in preparation (theory); 
we are grateful to P. Bertet for informing us about this work 



4 Periodically modulated quantum nonlinear oscillators 



drive an oscillator into the amplitude range where Auj > F. This makes underdamped 
oscillators advantageous for studying quantum phenomena far from equilibrium. 

7.2.1 Oscillator Hamiltonian in the rotating frame 

The most frequently used types of resonant modulation of an oscillator are modulation 
by a resonant additive force A cos upt with frequency cop close to wq and parametric 
modulation by force Fcoscjpt with cop close 2wo- For moderately strong modulation 
it is often sufficient to take into account only the leading-order oscillator nonlincarity 
which leads to the amplitude dependence of the vibration frequency. It is sometimes 
called Kerr nonlincarity, and the corresponding model of the oscillator is called the 
Dufhng model. The Hamiltonian of the Duffing oscillator is 

Ho = ^p' + ^cjiq' + ^^q* + Hp{t), (7.1) 

where q and p are the oscillator coordinate and momentum, the mass is set equal to 
one, and 7 is the anharmonicity parameter. 

The modulation term Hp in eq. (j7.ip for additive {Hp ~ -ffadd) and parametric 
(Hp = i/par) modulation has the form 

Hi,dd ^ -qAcosujpt, Hp^^ = ^q'^F cos ojpt. (7.2) 
The conditions that the modulation is resonant and not too strong are 

|(5w|<wo, (5w = a;M-wo; l7l(9^)<'^o- (7-3) 

Here, wm is equal to ujp and LUp/2 for additive and parametric modulation, respec- 
tively; this is the frequency close to the oscillator eigenfrequency for these types of 
resonant modulation. For concreteness, we assume 7, > 0; for additive driving, the 
oscillator can be bistable for ^Suj > 0; we assume 6uj > for such driving. 

It is convenient to change to the rotating frame using the standard canonical trans- 
formation U{t) = exp (— ia^awM^), where a' and a are the raising and lowering oper- 
ators of the oscillator. We introduce slowly varying in time dimcnsionless coordinate 
Q and momentum P, using as a scaling factor the characteristic amplitude of forced 
vibrations C, see Table [TTI 

U\t)qU{t) = C{Q cosip + Pshlip), 

U\t)pU{t) = -Ca;M(Qsin^ - Pcos^). (7.4) 

For additive and parametric modulation Lp = </3add = ujpt and ip = <Ppar = {i-upt + TT)/2^ 
respectively. The commutation relation between P and Q has the form 

[P,Q]^-iX, X = h/{iOMC^). (7.5) 

Parameter X cc H plays the role of the Planck constant in the quantum dynamics in the 
rotating frame. It is determined by the oscillator nonlincarity, A oc 7, see Table 17.11 
For characteristic |(5|,|P| < 1, where (g^) < C^, the last inequality in eqn (|7.3p 
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Table 7.1 Parameters of a resonantly modulated oscillator 



Additive driving 



Parametric driving 



Amplitude scale 


C = 




C = 




Scaled Planck constant 


A = 


ih'y /8lup{ll>f — ^o) 


A 


3fi^^i|7/F| 


Control parameter 


/?- 


37ylV32a;|,(a;F - wo)^ 


M = 


wf(wf - 2u!q)/\F\ 


Scaled decay rate* 


K = 


fl^^ = T/\luf — (^o\ 


K = 


= 2Tluf/\F\ 



* Notations Q and were used in some of our previous papers, cf. pykman (2007| l 

coincides with the first inequality in this equation for additive modulation, whereas 
for parametric modulation is gives condition <i; Wq. 

In the range (|7.3p the oscillator dynamics can be analyzed in the rotating wave 
approximation (RWA). The Hamiltonian in the rotating frame is 



Ho = U^HoU - ihU^U « {3E,i/8)g, E,, = 



(7.6) 



where Es\ ^ l{<l^) is the characteristic energy of motion in the rotating frame. This 
motion is slow on the time scale w^^. Operator g = g{Q, P) in eqn (|7.6p is independent 
of time. For additive and parametric modulation, respectively, we have 



and 



gadd{Q,P) 



p2 + 02_l)2_^l/2, 



l{P' + Q'T + l{i-i^)P'-lil + i.)Q' 



(7.7) 



(7.8) 



In Fig. 17.21 we show gadd{Q,P) and (?par(Q,-P) as functions of classical coordinate 





Fig. 7.2 (Color) The dimensionless Hamiltonian functions of the oscillator for additive (left 
panel) and parametric (right panel) modulation. The plots refer, respectively, to /3 = 0.01 
and to ^ = —0.1. In the presence of weak dissipation, the minimum and the local maximum 
of (7add and the minima of gipar become classically stable states of forced vibrations in the lab 
frame. 



and momentum. Each of these functions depends on one dimensionless parameter, 
f3 and /x, respectively, which characterizes the ratio of the modulation strength to 
the frequency detuning. These parameters are given in Table [TT] For < /3 < 4/27 
function (/add has the form of a tilted Mexican hat, with a local maximum and with 
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a minimum at the lowest point of the rim. In the presence of weak dissipation these 
extrema correspond to classically stable states of forced vibrations with small and 
large amplitude, respectively, see Fig. IT-lf a). 

For — 1 < ^ < 1, function ^par has two minima, which in the presence of weak 
dissipation correspond to stable vibrational states, see Fig. I7.1f b). Function gpar has 
symmetry gpar((5, -P) = .9pai (~Q, ~-P)- This is a consequence of the time-translation 
symmetry Ho{t) = Ho{t + 2tt/ujf), as seen from eqn (|7.2p and (|7.4p . Respectively, the 
vibration amplitudes in the stable states are the same, but the vibration phases differ 
by TT, characteristic of parametric resonance. 

7.2.2 Quasienergy spectrum 

Operator g plays the role of dimensionless Hamiltonian of the modulated oscillator in 
the rotating frame. In the RWA, the Schrodinger equation in dimensionless slow time 
T reads 

iA?A EE iXdrip ^gil^, r = tX-fC^/h = {XEsi/h)t; (7.9) 

T = t5uj and r = tF/2ujp for additive and parametric modulation, respectively. 

Operator g has a discrete spectrum, g\n) = gn\n). The eigenvalues gn have sim- 
ple physical meaning. A periodically modulated oscillator does not have stationary 
states with conserved energy in the lab frame. It is rather described by the Floquet, 
or quasienergy states = U{t)ipe{T), which is a consequence of the periodic- 

ity of the Hamiltonian HQ{t) = Ho{t + 2tt/ujp). One can seek a solution of the full 
Schrodinger equation ihdt'^ ~ Ho{t)'i' in the form + tu) = cxp(— ietM/^)^'e(0: 
where = Ztt/wm (^m = 2t: /lof and = 4:t:/lof for additive and parametric 
modulation, respectively). This expression defines quasienergy e. We note that, for 
parametric modulation, we use the doubled modulation period when defining e; this 
is convenient for the description of period-two states of the oscillator. 

From eqns ()7.2p and ()7.9p . in the RWA the oscillator quasienergies are simply 
related to the eigenvalues of g, £„ = {3Esi/8)gn, i.e., gn is a scaled quasienergy. Here 
we are using an extended e-axis rather than limiting e to the analog of the first Brillouin 
zone < e < Hujm- The scaled quasienergy spectrum for parametric modulation in the 
neglect of tunneling is sketched in the right panel of Fig. 17.31 In the region of bistability, 
— l<fi<l, the states with g„ < are degenerate. Function gpmiQjP = 0) has a 
form of a symmetric double- well potential, and the states with gn < remind intrawell 
states of a particle in such potential. Semiclassically, in the rotating frame the oscillator 
moves along orbits, which are cross-sections of the surface gpai in Fig. 17.21 by planes 

The structure of quasienergy states of an additively driven oscillator can be un- 
derstood in a similar way by thinking of the cross-sections of the surface g^ddiQ, P) 
in Fig. 17.21 bv planes (?add = dn- The eigenstates localized near the local maximum 
of ga,dd{Q,P) correspond to semiclassical orbits on the surface of the "inner dome" 
of 5add(Q, P)', these states become stronger localized as g„ increases toward the local 
maximum of g^ddiQ, P)- This is in contrast with the conventional picture of a particle 
in a potential well, where the localization becomes stronger with decreasing energy. 

Many interesting and unusual features of the dynamics described by eqn ()7.9|) 
follow from the fact that functions (7add : ffpar do not have the form of a sum of the 
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Fig. 7.3 Panel (a); oscillator relaxation and excitation. Relaxation is due to transitions 
between the Fock states with energies En ~ h!jJo{N + 1/2) accompanied by emission (or 
absorption, for nonzero temperature) of excitations in the bath, e.g., photons. The station- 
ary vibrational state is formed on balance between relaxation and excitation by periodic 
modulation F{t). Panel (b): the effective RWA Hamiltonian of the parametrically modulated 
oscillator .gpar(Q,f = 0) for fj, = —0.1, with sketched quasienergy levels. The dashed arrows 
indicate transitions between quasienergy states, which are due to emission of excitations in 
the bath. The solid arrows indicate the change of g in relaxation and quantum heating. 

kinetic and potential energies. These features are seen, in particular, in tunneling, 
which is significantly modified compared to the conventional picture, because for 
a given g{Q, P) the momentum P as function of coordinate Q has 4 rather than 
2 branches ( [Dmitriev and Dyakonov, 19866 Serban and Wilhelm, 20071. One of the 



consequences is that, for example, for a parametric oscillator, decay of the wave func- 
tion in the classically inaccessible region of Q can be accompanied by oscillations, 
leading to undcr-barrier interference ( [Marthaler and Dykman, 2007[ ). 

As seen from cqn (j7.9p . the quasienergy level spacing is oc XEs\. From cqn (|7.3|) . it 
is small compared to the distance between the oscillator energy levels in the absence 
of modulation, |e„ — ^ ^Esi <^ huo- 

7.2.3 Qualitative picture of quantum heating and the master equation 

Quantum heating is most easy to understand in the case where the oscillator decay 
rate F is small not only compared to wq, but also to the distance between quasienergy 
levels in frequency units, F <^ XEsi/h. In this Section we consider relaxation due to 
coupling to a bosonic thermal bath with the coupling Hamiltonian Hi linear in the 
oscillator coordinate and momentum and thus in the oscillator ladder operators a, 



(Schwinger, 19611 



H, = ahb + U.c, r = T{ujQ) = h-^Rc / dt{[hl{t), hi,{0)])be"^°\ (7.10) 





where hh depends on the bath variables only and (. . .)b denotes thermal averaging 
over the bath states; in what follows we assume (/ib)b = 0. 

For the interaction (|7.10p . in the absence of modulation relaxation is due to tran- 
sitions between adjacent oscillator Fock states \N). For zero bath temperature, these 
transitions occur only downward in energy, with emission of excitations in the bath, see 
the left panel of Fig. 17.31 For a smooth density of states of the bath, resonant modula- 
tion does not change the decay rate, F(wm) ~ r(a;o)- However, it excites the oscillator. 
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as sketched in Fig. 17.31 In the stationary vibrational state the energy provided by the 
modulation is balanced by relaxation. 

As outlined in Sec. 17.11 the randomness of the transitions, i.e., the quantum noise 
that accompanies relaxation, leads to quantum heating. The heating can be understood 
from Fig. 17.31 bv noticing that the quasienergy states \n) {g\n) = gnln)) sketched in 
the right panel are linear combinations of the Fock states |A^) in the left panel, |n) = 
X^AT '^"A'l-^)- Therefore transitions between the Fock states downward in oscillator 
energy correspond to transitions both downward and upward in quasienergy, with 
different rates. If the state with minimal g„ is the stable state, transitions downward are 
more likely, but upward transitions still have nonzero rates. The outcome is diffusion 
over quasienergy states away from the minimum of g (or the extremum of g, for 
an additively modulated oscillator, see Fig. 17. 2[) . that accompanies drift (relaxation) 
toward the minimum (extremum) of g. 

For a thermal equilibrium system with nondegenerate energy levels Ea, the ratio 
of the rates of interstate transitions \a) — >■ |/3) and \(3) — ^ \a) due to weak coupling to 
a bath is W^^/W^ = exp[(i^Q, — Ep)/kBT] and is fully determined by temperature. 
Similarly, the ratio Wnm/Wmn of the transition rates between quasienergy states \n) ^• 
I to) and I to) — >■ |n) characterizes the effective temperature Te of the distribution over 
these states. It is nonzero even where the bath temperature T = 0. 

A modulated oscillator does not have detailed balance for T > 0. The transitions 
\n) — >■ I to) are not limited to to = n ± 1 (in which case detailed balance would hold 
automatically), and the ratio Wnm/Wmn cannot be written as exp[(e„ —£m)/Te\ with 
the same effective temperature Te for all n and to. In other words, the stationary 
distribution is generally not of the Boltzmann form, it can be described by an e- or, 
equivalently, ^-dependent temperature. 

A complete analysis of the distribution can be done using the master equation for 
the oscillator density matrix p. In slow dimensionless time r, for the coupling to a 
thermal reservoir (|7.10p this equation reads 

p=drP= iX^^lp^g] — kp^ kp = K,{h + l){a^ ap — 2apa^ + pa' a) 

+Kfi{aa} p ~ 2a) pa + paa}), k = hV / \Es\. (7-11) 

Here, the term cx [p, g\ describes dissipation- free motion, cf . eqn. ()7.9p . Operator kp 
describes dissipation and has the same form as in the absence of oscillator modulation 
( [Mandel and Wolf, 1995| |Dykman and Krivoglaz, 1984[ ), k is the dimensionless decay 
rate, see Table [7?T1 a is the lowering operator, and n is the oscillator Planck number, 

a= (2A)-^/2(Q + iP), n = n(wo) [exp(fta;o/fcsr) - (7.12) 

[in the lab frame, a has an extra factor exjp[—iujy[t)]] the oscillator eigenfrequency ujq 
is defined so as to incorporate the renormalization due to the coupling (j7.10p . 

7.2.4 Semiclassical dynamics condition 

Of utmost interest is the parameter range where quantum fluctuations about the states 
of the forced vibrations of the oscillator are small compared to the interstate distance 
in phase space and the oscillator motion in the rotating frame is semiclassical. In 
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terms of Fig. 17.31 it means that the number of states m the wells of gpai {Q, P) is large. 
Similarly, a large number of states are localized near the extrema of (7add('5, -P), for 
additive modulation. As seen from eqns (j7.7|) - (|7.9p . this requires the dimcnsionless 
Planck constant to be small, 

A<1. (7.13) 

In the range (|7.13p the rate of switching between the stable states of forced vibrations 
Wgw is exponentially small, — InWtun oc for T = 0, see Sec. 17.31 We will assume 
that VFsw ^ T- Then over time ~ F"^ the oscillator will reach the state of forced 
vibrations in the vicinity of which it was prepared initially and will then fluctuate 
about it; switching between the states takes an exponentially longer time, see Sec. 17.31 

7.2.5 Quantum temperature 

Near the extrema function g{Q, P) is parabolic. The oscillator motion in the rotat- 
ing frame is mostly harmonic vibrations about these extrema provided the quantum 
smearing cx X^/"^ is small compared to the scale in phase space where the nonparabol- 
icity of g{Q: P) becomes substantial (see also Sec. l7.4p . If we disregard dissipation, the 
motion is described by the Heisenberg equations 

Q = dQ/dT ^-i\-~\Q,g\, P = dP/dT = -iX-^[P,g] (F ^ 0) 

linearized in Q—Qo, P. Here, Qq, P = is the position of the considered extremum of g. 

1/2 

The dimcnsionless vibration frequency is vq — \gQQgpp\ [the derivatives of g{Q, P) 
are calculated for Q — Qo,P — 0]. Its dependence on the oscillator parameters is 
shown in Fig. 17.91 In the ground vibrational state the average values of {Q — Qo)^ and 
P^ are different, which indicates that forced vibrations in the lab frame are squeezed. 

It is convenient to change from Q — Qq , P to the appropriate raising and lowering 
operators and b using the standard squeezing transformation 

Q-Qo + iP= (2A)^/^(fecosh¥J* - 6^ sinh(^,), 

g^g{Qo,0) + Xi^o{b^b + l/2)sgngQQ, i^o ^ \gQQgpp\^^\ (7.14) 

which makes the mapping on small-amplitude vibrations explicit; in eqn (|7.14|) . tanh ip^ ~ 

i\9QQ\'/' - \9Pp\^/')/i\9QQ\'/' + \9PP\'/')- 

The distribution over quasicnergy states near the extrema of g{Q, P) can be found 
from eqn (j7.1ip . For small damping, k ^ z/q, off-diagonal matrix elements of p in the 
basis of quasienergy wave functions are small, \pnm\ <C PnmPmm for m ^ n. To find 
the diagonal matrix elements one should substitute into eqn (j7.1ip g from (|7.14p and, 
using eqns. ()7.12|) and ()7.14p . express a, a) in kp in terms of operators 6,6^ keeping 
only bilinear terms that contain both b and W while disregarding terms with 6^, (6^)^ 
as well as the terms linear in 6, 6^. Then operator kp in terms of 6, b^ becomes of the 
same form as in terms of a,a\ except that n is replaced with fie, 

fie = + (2n + 1) sinh'^ 7^. = Ai^o/ ln[(^e + l)/^e]- (7-15) 

The stationary solution of the resulting equation is of the Boltzmann type, p^'^'' = (?T-e + 
1)^^ eyiY){—XvQb^b/Te) ( |Dykman et ai, 2011[ ). Parameter Te given in eqn. (|7.15p is the 
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Additive driving, small-amplitude state 
n= 0.6 



Additive driving, large-amplitude state 




Parametric modulation 









0.3 








Fig. 7.4 The effective Plancli number fie, for vibrations about tlie state of forced vibra- 
tions. The left, central, and right panels refer, respectively, to the small- and large-amplitude 
states of an additively modulated oscillator and to a parametrically modulated oscillator. 
The damping is assumed to be small, k -^vc^va eqn (|7.11[) . 

efTective dimensionless temperature of vibrations about the stable state in the rotating 
frame. For fi = the result coincides vifith that for a driven oscillator resonantly coupled 
to a two- level system (Peano and Thorwart I2010ffl[ I20106p . The distributions over 
quasienergy states for other systems and other relaxation mechanisms were discussed 



by Verso and Ankerhold (20101 and Ketzmerick and Wustmann (2010). 

From eqn (j7.15| ). fif = sinh ip^, > even where the Planck number of the original 
oscillator is zero. In the opposite limit of high temperature, ksT ^ fn^JQ, we have 
Te <x T, and p'^'^ looks like a Boltzmann distribution of an oscillator with frequency 
oJo/{i + 2 sinh^ (f^) < ujo- The dependence of He on the parameter that characterizes 
a resonantly modulated oscillator in the small damping limit is shown in Fig. 17.41 

Expression (|7.15p is simplified also near the bifurcation point (the bifurcation pa- 
rameter value) where the corresponding stable vibrational state disappears. If dissipa- 
tion is disregarded, at the bifurcation point the corresponding extremum (two extrema, 
for parametric oscillator) and the saddle point oi g{Q, P) merge, whereas for the large- 
amplitude state for additive modulation, the values of ^(Q, P) at the extremum and 
the saddle point coincide. As the parameters /3 oi fi approach their bifurcation values 
Pb or fiB, see Fig. 17. 1[ ggg — > or, for the large-amplitude state of an additively 
modulated oscillator {(3 — )• /3bi), gpp — )• 0. Near a bifurcation point z/q ^ 1- Then 
Te ~ (A/4)(2f;, + 1)|5pp|; for the large-amplitude state of an additively modulated 
oscillator gpp should be replaced with ggg. As we see, Te displays a characteristic 
temperature dependence described by the factor 2n -I- I. To the leading order, it is 
independent of the distance to the bifurcation point ri = j3 — /Sb or rj — fi — hb- In 
contrast, fie oc sharply increases with decreasing This is because the system 
becomes "soft" near a bifurcation point, and respectively, the distribution broadens. 



The scaling of Ue with 77 is rie oc j/q ^ 



cx |r/| 



with — 1/2 for parametric modulation 



and = 1/4 for additive modulation. 

The results for Te apply for not too small |?7|. One constraint on \r}\ is classical, 
vq 3> The other is quantum: the dimensionless spacing of quasienergy levels Xvq cx 
\r]\^'^ should be small compared to the difference between the values of g at the 
extremum and the saddle point of g{Q, P) that merge at the bifurcation point; for the 
corresponding extremum, for parametric modulation Ag oc 77^, whereas for additive 
driving Ap oc rf/'^ for j3 — > /3s2 and Ag oc 77-^/^ for /3 — > Pbi- 



Quantum activation 11 



7.3 Quantum activation 

The diffusion over quasienergy states, tfiat underlies quantum lieating, also populates 
states near the top of the quasienergy barrier in Fig. 17.3( b). As a result, if the sys- 
tem was initially occupying the left well of g{Q, P), for example, it will diffuse to the 
barrier top and switch to the right well. Such an overbarrier transition reminds the 
conventional switching via thermal activation in systems close to thermal equilibrium 
( [Kramers, 1940p , where the states near the barrier top are populated as a result of 
thermal fluctuations. In the case of a modulated oscillator, for low temperature the 
effect is due to quantum fluctuations. Respectively, the switching mechanism can be 
called quantum activation. It applies to both parametrically and additively modu- 
lated oscillators: in the latter case switching occurs, with probability ^ 1/2, once the 
oscillator located initially near an extremum of g{Q, P) reaches the saddle, see Fig. 17.21 
For small effective Planck constant A, the switching rate is exponentially small, 

ocexp(-i?^/A). (7.16) 

This estimate can be easily understood from Fig. 17.3( b). If the ratio of the typical rates 
of transitions up and down in quasienergy is W^/W^ < 1 and the transitions occur 
primarily between a few neighboring quasienergy states, the population of the levels 
close to the barrier top is ^ (W^/Wi)^^ , where M is the number of intrawell states. 
Since the dimcnsionless level spacing gn+i—gn ^ ^ the well depth .gmax~5min ^ Ij 
we have M ^ A~^, which immediately gives — InWsw ~ A~^ for low temperatures. 
This estimate applies to both additively and parametrically modulated oscillators. 

Even though the effective activation energy for a given stable state Ra is deter- 
mined by quantum fluctuations for low T, quantum activation is not tunneling, it is the 
result of coupling of the oscillator to a thermal bath and the quantum noise that accom- 
panies relaxation. Therefore finding Ra requires solving the master equation (|7.11|) . For 
small A this can be done using the WKB approximation. The problem is qualitatively 
different from that of switching in systems close to thermal equilibrium, which can be 



efficiently approached using the instanton technique ( [Langer, 1967[ Coleman, 1977 



[Afffeck, 1981 Caldeira and Leggett, 19831. A modulated oscillator does not have de- 



tailed balance, generally, and its distribution is not of the Boltzmann form and is not 
characterized by the partition function. 

For small relaxation rate, switching of a modulated oscillator can also occur via 
dynamical tunneling with constant quasienergy, for example via tunneling between 
equal-quasienergy states in the left and right wells of g{Q,P) in Fig. I7.3f b). see 
Larsen and Bloemb crgcn 1976':'Sazonov and Finkclst cin 19761 [Dmit ricv a-pd Dyakonov 19866* 



Vogel and Risken 1988; ,Wiclinga and Milburn 1993j Peano and Thor wart 2006: ,Marthalcr and Dyk 



Serban and Wilhelm 20071 As we explain below, switching via tunneling becomes sub- 



stantial only where the relaxation rate of the oscillator is exponentially small, otherwise 
switching occurs via quantum activation. 

7.3.1 The WKB svi^itching paths 

The rate of switching from a given vibrational state can be found from the quasi- 
stationary solution of the master equation. It is formed on times F^^ ^ t <C 
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in a broad range of phase space. As mentioned previously, in this time domain, the 
oscillator prepared at i = near the considered stable state in phase space will have 
come to local equilibrium, but most likely will not have switched to another state. 

The physical picture of switching can be understood from the classical phase por- 
trait for the range of bistability shown in Fig. 17.51 The positions of the stable states 
on the oscillator phase plane in the rotating frame {Qa.,Pa.) are given by the stable 
solutions of classical equations of motion 

Q = dpgiQ, P) ^kQ, P = -5q5(Q, P) - nP (7.17) 

These equations immediately follow from master equation (|7.11|) if in the equations for 
the average values (Q) = Tr Qp, (P) = Tr Pp one disregards fluctuations, {Q^^P"^) ^■ 

(Q)"(p)™ ^ g"P". 




Fig. 7.5 The phase portraits of the classical additively (a) and parametrically (b) modulated 
oscillators in the rotating frame. The attractors ai,2 correspond to the stable states of forced 
vibrations in the lab frame. Their positions (Qajfa) give the scaled vibration amplitudes 
^osc = (Qa + PaV^^ in Fig. 17.11 The separatrix that goes through the saddle point S is the 
boundary of the basins of attraction to different attractors. The phase portrait in (b) has 
inversion symmetry. The parameters in (a) are /? = 1/27, k — 0.15; in (b) p, = 0.2, k — 0.3. 

For t <C the oscillator is mostly localized in phase space near the initially 

occupied stable state, with localization length ^ A^^^ for low temperatures. The tail of 
its quasi-stationary distribution is formed by large rare fluctuations. In the semiclassi- 
cal picture, switching occurs as a result of a large fluctuation that brings the oscillator 
to the basin of attraction of the initially empty stable state. From there, the oscillator 
will most likely move to this state staying within oc A^/^ from the classical trajectory 
(|7.17p . The rate W^sw is determined by the maximal rate of an appropriate fluctuation. 

In this section we will be interested in finding the switching exponent, i.e., the 
leading-order term in In Wsw ■ This can be done ( [Dykman and Smelyansky, 1988[ ) us- 
ing the density matrix in the coordinate representation p{Qi, Q2) = {Qi\p\Q2)- Equa- 
tion (|7.1ip in this representation has the form 

p{Qi,Q2) = -i>^^^'H{Qi,Q2,-i\dQ^,-i\dQ^)p{Qi,Q2), 

where 



n{Qi,Q2, Pi, P2) = giQuPi) - 9{Q2, P2) - k{PiQ2 + P2Q1 - iX) 

-^iKi2n + 1) [(Qi - Q2)' + {Pi + . (7.18) 
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Along with p{Qi, Q2) it is convenient to consider the density matrix in the Wigner 
representation 

PwiQ, P) = j d^e-'^P'^p{Q + Q - i/2). (7.19) 

For the oscillator, function pw in the quasi-stationary regime has a Gaussian peak at 
the initially occupied stable state (Qaj -Pa); this form of pw is generic for semiclassical 
systems (for k ^ z/q it follows from the results of Sec. I7.2.5p . It rapidly falls off away 
from Qa, Pa- The switching exponent is determined by the maximal value of pw{Q, P) 
for Q, P inside the basin of attraction of the initially empty stable state or on the 
basin boundary, i.e., it is determined by the tail of pwiQ^P) a-nd in turn, by the tail 
of p{Qi,Q2)- 

In the spirit of the WKB approximation, one can seek p{Qi, Q2) on the tail in the 
eikonal form. To the leading order in A in the quasi-stationary regime {p = 0) wc have 

P(gi,g2) =exp[i5(Qi,g2)/A], n{QuQ2,dQ,S,dQ,S)^Q. (7.20) 

Equation (|7.20p has the form of the Hamilton- Jacobi equation for an auxiliary classical 
system with coordinates Qi, Q2 and action S, with equations of motion 

Qj=dp^HiQi,Q2,Pi,P2), Pj=-dQ^H{Qi,Q2,PuP2) (j = l,2). (7.21) 

Equations (fr20| and map the problem of finding the tail of the density 

matrix of a dissipative quantum oscillator onto the problem of finding classical Hamil- 
tonian trajectories. The initial conditions for these trajectories follow from the Gaus- 
sian form of Pw near Qa,Pa- From ()7.19p and ()7.20p . for the trajectories coming 
from the stable state of the oscillator the initial conditions are Qi ~ Q2 ^ Qa and 
Pi = —P2 = Pa, where Pi. 2 = dq-^ ^S. The analysis shows that the final point on the 
switching trajectory (|7.21|) is Qi ~ Q2 — Qs and Pi = — P2 = Ps, where (^5,^5) is 
the saddle point of the oscillator. 

The effective activation energy for switching from a given stable state is determined 
by action S calculated along the switching trajectory. 



Ra = Ini 



/oo 
-00 ' 



We took into account that the initial and final points on the trajectory are stationary 
states of the auxiliary system, therefore the time integral goes from —00 to 00. 

The formulation (|7.18p - (|7.22p reminds the conventional instanton formulation. 
The major distinction is that the motion occurs in real rather than imaginary time. At 
the same time, the switching trajectory is complex, as are also H and S. It is instructive 
to compare this formulation with the theory of noise-induced switching in classical 
systems. There, for Gaussian noise, the switching rate displays activation dependence 
on the noise intensity {ksT, for thermal noise). The effective activation energy can be 



calculated as action of an auxiliary Hamiltonian system (Freidlin and Wentzell, 1998 
|Dykman, 1990[ [Kamenev, 2011[ ). The trajectory followed by the auxiliary system is 
real and gives the most probable trajectory that the initial noise-driven system follows 
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in switching. Such trajectory has been seen in experiment, see the chapter by H. B. 
Chan and C. Stambaugh in this book and references therein. 

One can see that, due to the symmetry of the Hamihonian H, eqn (|7.2ip has 
a solution Q2{t) = Ql{t), P2{t) = —Pl{t), which satisfies the boundary conditions. 
This solution was used to find the activation energy Ra for an additively modulated 
oscillator as function of /3 for small damping ( [Dykman and Smelyansky, 1988[ ). 

7.3.2 Balance equation approach 

We now discuss an alternative approach, which immediately gives the switching rate 
and the distribution of the oscillator in the small damping limit. If the broadening of 
the quasienergy levels is small compared to the interlevel distance, one can disregard 
off-diagonal matrix elements pnm {n 7^ hi the basis of quasienergy wave functions. 
Then eqn (|7.1ip is reduced to a balance equation for state populations 

P«n = {WynnPrnm ~ WnynPnn) , 1V™„ = 2k [(fl + 1) | a„„ ^ + n| a™„ H , (7.23) 

where a„,„ = (n|a|m) (we remind that a is the oscillator lowering operator). We 
disregard tunneling when defining functions |n) = '0„((5), i.e., we use the "intrawell" 
wave functions in Fig. I7.3r b): the effect of tunneling is exponentially small for A <C 1. 

Matrix elements a,„„ can be calculated in an explicit form using the WKB ap- 
proximation. It relates the problem to that of classical conservative motion Q = 
dpg,P ~ —dqg. Such motion is periodic oscillations in time Q{t; g)^ P{t; g) with 
given g{Q, P) = g and with dimcnsionlcss frequency v{g) that depends on g. For not 
too large |m — n|, matrix clement is given by the (m — ri)th Fourier component of 
the periodic function a(r;(7„) = (2A)~^/^ [Q{'^'ig) +*-?'(■'"; ff)] calculated for the classi- 
cal orbit g{Q, P) — .g„. Formally, we require that |m — n| ^ n, but the results apply 
also near the cxtrcma of g{Q,P) where n ^ 1, since ipniQ) are close to the wave 
functions of a harmonic oscillator for such n. 

The evaluation of amn simplifies if one notices that g{Q, P) is quartic in Q, P and 
even in P. Because of that, the orbits Q{t; g), P{t; g) are described by the Jacobi el- 
liptic functions and are double periodic in t. To calculate one can then integrate 
a(r; gn) exp[iT(n — m)z^((7„)] along an appropriately chosen closed contour on the com- 
plex r-plane. The result is determined by the pole of a(r; gn) and has a simple form 
( Marthaler and Dykman, 2006| ). In particular, Wmn exponentially decays with |m — n| 
for 1 <C |to — n|, with the exponent that depends on the sign of to — n. 

The quasi-stationary populations of neighboring states n and n ± 1 in eqn (|7.23p 
can differ significantly. However, as we will see, lnp„„ is a smooth function of n. 
Respectively, we seek the quasi-stationary distribution in the eikonal form, pnn = 
exp(— i?,i/A), Rn = R[gn)- To the leading order in A, eqn (|7.23p then reads 

Wn+kn {1 - exp [-ky{gn)R'{gn)]} = 0, R'{g) = dR/dg, (7.24) 

where we used gn+k ~ 5« + \kiy{gn), Rn+k ~ Rn + Xkiy{gn)R'ign) and Wnn-k « 
Wn+kn for |fc | < n 

From eqn (|7.24|) . R'{g) is independent of A and is given by a solution of a polynomial 
equation. For g close to its value gs, = giQa, Pa.) at a stable state (Pa — )■ for k — >• 0), 
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the solution of eqn (|7.24p is of the Boltzmann form, p„„ cx exp[— At'o?ii?'(5a)]; as one 
can show, it coincides with the result of Sec. 17.2.51 with R'{ga) = However, this is 

only the asymptotic solution, generally the distribution is not described by an effective 
temperature, because R'{g) varies with g. We note that the corrections disregarded in 
deriving cqn (|7.24[) arc cx A, which justifies this equation for A ^ 1. 

The effective activation energy of switching from a given stable state is determined 
by the quasi-stationary occupation of states near the saddle point, with dimensionless 
quasicnergy gs = g{Qs, Ps)- Therefore, to the leading order in A 



Ra 



dgR\g). 



(7.25) 



In Fig. 17.61 we show the activation energy of switching between period-two states of a 
parametric oscillator obtained from eqns (|7.24p and (|7.25p ( [Marthaler and Dykman, 2006p . 
For small damping, Ra depends on two parameters, the scaled frequency detuning of 
the modulating field fx and the Planck number of the oscillator n. As seen in the right 
panel, the value of (2n -I- \)Ra decreases with increasing temperature and already 
for n = 1 becomes very close to the result for the classical range fi ^ 1. In this 
range, switching is thermally activated an d Ra (x 1/T JPykman et ai, 19981 ); cf- the 
experiments by Lapidus et al. (1999) and Chan and Stambaugh (2007). Both in the 
classical limit and near bifurcation points \i'(a)R'(a)\ «: 1 and eqn ( j7.23p becomes a 
linear equation for R'{g). 





Fig. 7.6 The tunneling exponent S'tun for tunneling between the extrema of g{Q,P) of a 
parametric oscillator and the switching activation energy Ra in the limit of small damping 
for different values of the Planck number ( Marthaler and Dykman, 2006[ ). The cross indicates 



the value of 5tun obtained by Wielinga and Milburn (1993 I. 



Near the bifurcation point in Fig. 17.11 where the stable period two states 
merge together and disappear, Ra {2h + — ^bi)^- For additive driving, for 

switching from the large and small-amplitude vibrational states in Fig. 17.11 near the 
corresponding bifurcation points Ra oc (2fl-|-l)~^(/3— and Ra oc (2n-|-l)^^(/3_B2 — 
^•j3/2^ respectively ( [Dykman and Smelyansky, 1988[ ). This scaling applies not too close 
to bifurcation points, where still z/q = i^iga.) 3> n. 



Detailed balance for T = 0. An important feature of the scmiclassical matrix el- 
ements, which follows from the double-periodicity of a{T;g), is that \amn/an 



^nm I ' 
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exp[(m — n)c{gn)] for |m — n| ^ n, with c{gn) that smoothly depends on n. There- 
fore for T = the system has detailed balance: the ratio of the transition rates is 
path- independent, WnrnW^k / {WkmWmn) = Wnk/Wkn, and eqn (|7.23p has a quasi- 
stationary solution Pmm/ Pnn = Wnm/Wmn- This gives the quasi-stationary popula- 
tion of states near the saddle point {Qs,Ps) relative to that near the initially occupied 
stable state (Qa, Pa), and thus the activation energy Ra = — J^^ dgv~^ {g)c{g) . For a 
parametric oscillator Ra for T = is shown in Fig. 17.61 

A remarkable property of the detailed balance solution seen from Fig. l7.6l is fragility: 
the value of Ra for T = differs from Ra for T 0. The fragility emerges for small 
A in the limit where the relaxation rate F oc k — ^ 0. For nonzero k, the transition from 
T = to nonzero T solutions should be continuous ( [Dykman and Smelyansky, 1988 ); 
more work is required to study this transition. The detailed balance condition for T = 
applies for arbitrary k. This was used to find the stationary probability distribution for 
both additively and parametrically modulated oscillators ( [Drummond and Walls, 1980[ 



Kryuchkyan and Kheruntsyan, 1996) and the switching rate in the overdampcd limit 



for parametric modulation (Drummond and Kinsler, 1989). 

We now compare switching via quantum activation and dynamical tunneling. The 
rate of switching via tunneling is oc exp(— 2S'tun/A), where S'tun is the tunneling ac- 
tion near the corresponding extremum of g{Q,P). The tunneling rate prefactor is 
(X voXEsi/h. The prefactor in the rate of switching via quantum activation for small 
damping is oc F oc k. From Fig. 17.61 and from similar results for additively modulated 
oscillator, 25'tun > Ra ( |Dykman and Smelyansky, 1988|[Marthaler and Dykman, 2006] ). 
Therefore unless k is exponentially small, oscillator switches via quantum activation 
rather than tunneling. 



7.3.3 Switching rates near bifurcation points 

The analysis of switching near bifurcation parameter values is particularly impor- 
tant. In this range the switching rates display universal, model-independent features; 
also, the range is interesting for many applications, in particular, for the Joscphson 



bifurcation amplifiers (Vijay et ai, 2009). Since the frequency of vibrations about a 
mctastable state i>q rapidly decreases with the decreasing distance to a bifurcation 
point ?7, see Sec. 17.2.51 for small ry the oscillator motion in the rotating frame is often 
overdamped, i/q k, which we will assume in this subsection to be the case. 

Near a bifurcation point, the behavior of the system is controlled by a "soft mode" , 
a dynamical variable that slowly changes in time. Without fluctuations the occurrence 



of such variable is well-known in classical dynamics ( Guckenheimer and Holmes, 1997 ) . 
It emerges because, for small rj, the stable and unstable states of the system are close 
to each other in phase space. In Fig. I7.5f a) the attractor corresponding to stable 
forced vibrations with large or small amplitude becomes close to the saddle point (the 
saddle-node bifurcation), whereas in Fig. l7.5f b) for small ^—ubi the attractors become 
close to each other and to the saddle point between them (the supercritical pitchfork 
bifurcation). In appropriately scaled variables, slow motion along the direction between 
the close states is described by equation x = —U'{x). The potentials U{x) for the 
bifurcations of interest are sketched in Fig. 17.71 
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Fig. 7.7 EfTective potentials for overdamped motion near (a) tiie saddle-node and (b) the su- 
percritical pitchfork (onset of stable period-two vibrations) bifurcations. The equation of mo- 
tion in scaled variables is i = —U'{x) with U{x) = —x^/3+rix in (a) and U{x) = x* /4 — rix^ /2 
in (b); the plots refer to 7; = 1. The minimum and maximum of U in (a) correspond to stable 
and unstable states of forced vibrations for additive modulation; the minima in (b) corre- 
spond to the stable vibrations with opposite phase, for parametric modulation. The specific 
form of U{x) is obtained by keeping the lowest order terms in x and rj compatible with the 
condition of merging of one or two stable and an unstable state of the system for 77 = 0. 

A theory of the switching rate based on the master equation in the Wigner repre- 
sentation was developed earlier ( [Dykman, 2007 1. Here we sketch a somewhat simpler 
derivation based on the quantum Langevin equation. Even without writing this equa- 
tion one can sec that, since the motion is overdamped and is characterized by one slow 
dynamical variable with no conjugate variable, commutation relations for this slow 
variable are irrelevant. Then its fluctuations are the same as in the case of a classical 
oscillator ( [Dykman and Krivoglaz, 1980 1, the only difference being that the fluctua- 
tion intensity is determined by quantum rather than classical noise from the thermal 
bath. From Fig. 17.31 for nonzero bath temperature the noise intensity is proportional 
to the overall rate of the bath-induced transitions up and down between the Fock 
states of the oscillator, which in turn is proportional to 2n-|- 1 (the Einstein relation). 
Respectively, using the quantum to classical correspondence for high temperature, in 
the expression for the rate of classical thermally activated switching one should replace 
kgT with (2n 4- l)fta;o/2, which indeed gives the right answer. 

The quantum Langevin equation in the rotating frame is an extension of the Heisen- 
berg equation for Q,P that includes the effect of coupling to a thermal bath. In the 
same approximation that led to the Markov master equation (j7.1ip . in slow time r 



g = -tX-'[Q,g] -kQ + /q(t), P = -t\-'[P,g] - kP + fp{T). 



(7.26) 



Here, fq.p are quantum noise operators. Equation (|7.26p is well known ( [Ford et al., 1965| 
for a harmonic oscillator linearly coupled to a bath of harmonic oscillators. It applies 
also for a more general form of the coupling to the bath Hi , eqn (|7.10p . To the leading 
order in Hi it can be obtained just by iterating the Heisenberg equations of motion 
for the bath ( |Lax, 1966| . The nonlinearity of our oscillator is relatively weak and does 
not affect the form of the dissipative and noise terms in eqn ()7.26|) . 

From eqn (|7.10[) . /q,/p are linear combinations of operators /ib(i) exp(— iwm*), 
hl^{t) exp{iu!Mt) calculated disregarding the coupling to the oscillator. For a smooth 
around wm power spectrum of /ib, the noise is (5-correlated in slow time. 



(/Q(r)/Q(r'))b = (/p(r)/p(T'))b = AK(2n+ l)<5(r - r'), 



(7.27) 
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and ([/q(t), /p(T')])b = 2i\K6{T — r'). This commutation condition guarantees that 
the commutation relation [Q, P] = iX does not change in time. The noise correlators 
are understood here in the Stratonovich sense ( [Van Kampen, 2007[ ); in particular, 
([/QM,P(r)])b = ([g(r),/p(r)])b=zA/c. _ 

Near a bifurcation point we can simplify eqn (j7.26p using essentially the same 
approach as for classical systems ( [Dykman and Krivoglaz, 1980 ). We can change to 



operators Q ~ Qb and P — Pb ; here Qb and Pb are the classical values of Q and P 
at the bifurcation point given by the appropriate stationary solutions of the classical 
noise- free equations of motion (|7.17|) . We can then further change from Q — Qb, P — Pb 
to Q' , P' by rotating coordinates in the [Q, P) plane. We choose the angle of rotation 
in such a way that at the bifurcation point the equation for P' does not contain linear 
in Q' , P' terms. The equation of motion for Q', on the other hand, has the form 

Q' = -AqqQ' + AqpP' + (nonlinear terms iuQ', P') + Jq,. 

Therefore, for small expectation values of Q',P', the operator P' is slowly varying 
in time compared to Q' . We note that, for an additively modulated oscillator, the 
rotation is not needed, P' = P — Pb, and for a paramctrically modulated oscillator 
Qb = Pb= 0. 

We will study the slow dynamics of P' for a small deviation rj of the control 
parameter from its bifurcation value. For \r]\ <^ 1, in dimensionless time 1/Aqq ^ 1 
operator Q' reaches its adiabatic form sa {Aqp/App)P' + 0{r]), while P' remains 
unchanged. The expression for Q' can be then substituted into equation for P' . The 
resulting equation for P' reads 

P' « Bov + BiT^P' + B^P'^ + B^P'^ + fp\T). (7.28) 

Here, near the pitchfork bifurcation point for a paramctrically modulated oscillator 
-Bo = -02 = by symmetry and i?i^3 ~ 1; near the saddle- node bifurcation point 
for an additively driven oscillator So,2 ~ 1 while the terms oc rjP\P'^ can be dis- 



regarded, cf. Guckenheimer and Holmes (1997). The explicit form of -Bo. 1,2. 3 follows 
from eqns 

The correlator of f'p is the same as of fp in eqn (|7.27p , whereas /g/ drops out from 
eqn (|7.28p . to the leading order in 77 {fgi enters the equation for P' with a coefficient 
(X P'). Therefore the 5-correlated noise fpi behaves as classical, as it commutes with 
itself, and then P' behaves as a classical variable. However, the intensity of fpi is 
(X 2n -|- 1, so that the fluctuations still have quantum origin. 

Equation (|7.28|) can be written as P' = —dp'U{P') + /p'(t). It maps the oscil- 
lator dynamics onto the dynamics of an overdamped classical Brownian particle in a 
potential U. The form of the potential depends on the nature of the bifurcation. For 
an additively modulated oscillator, where the stable and unstable states of forced vi- 
brations merge for (3 = /3b, this potential in rescaled variables is shown in Fig. I7.7r a). 
whereas for a paramctrically modulated oscillator where the period-two states merge 
for /X = /^si, it is shown in Fig. 17.7( b): the dynamics for small fi — fJ.B2 is described by 
a potential of the opposite sign. 
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Table 7.2 Quantum-activated switching near bifurcation points, W^sw = Slsw exp(— _Ra/A) 



Additive driving 



Parametric driving 



1 + 9k2=f (1-3k^ 



V 3/2 



Bifurcation points Pbi,2 = ^ 
Squared amplitude 

at bifurcation points (Aosc)bi,2 = ^ [2 ± (1 - 3k^)^/^] 

Distance to bifurcation rj ~ (3 — Pb 



Activation energy i?^ 
Prefactor flsw 
Auxiliary parameter 



2^2 



7^H3/V(2n+1) 



3k|6|i/2/3^ 

|<Sc.|(6,7/2)VV^/3V4 

b = p'J'[3iAljB~2]/2K^ 



i.2 = T(l-«^)^/^ 



(^osc)b1,2 = 
1] = II- HB 

\\^JiBW|{2n + l) 
r|r,UR|[l±Mi££)l 



Switching from a metastable vibrational state corresponds to a quantum-activated 
escape from the corresponding minimum of the potential U . The rate of escape via 
quantum tunneling is exponentially smaller ( [Dykman, 2007[ ) . The activation exponent 
Ra is determined by the height of the potential barrier, and therefore it displays a 
characteristic scaling dependence on the distance to the bifurcation point r/ = j3 — Pb 
or rj = fi — fj,B- The results are summarized in Table 17.21 

The scaling behavior of In Ws^ with the distance to the bifurcation point for clas- 
sical oscillators has been seen for additive driving near the saddle-node bifurcation 
points ( [Siddiqi et al., 2006} |Stambaugh and Chan, 2006 a[ ) and near the critical point 
(/3 = 8/27, K = 1/V3) where both stable vibrational states and the unstable state 
in Fig. FTSl a) merge together ( |Aldridge and Cleland, 2005[ 11, and for parametrically 
modulated oscillators near the pitchfork bifurcations ( [Chan and Stambaugh 2007 



see 

also the chapter by Chan and Stambaugh in this book). Recently the scaling behav- 
ior near the saddle-node bifurcation and the characteristic temperature dependence 
InVKsw oc (2n -|- were found also in the quantum regime ( jVijay et al. 2009} see 
also the chapter by R. Vijay ct aim this book). The results provide direct evidence in 
support of the mechanism of quantum activation. 



7.4 Power spectra of modulated quantum oscillators 

Of significant interest are spectra of a resonantly modulated oscillator, including the 

power spectrum and the spectrum of the response to an additional field ( [Dykman and Krivoglaz, 1979| 

[Drummond and Walls, 1980[ [CoUett and Walls, 1985[ ). Among other characteristics, 

the power spectrum determines the emission spectrum of the oscillator and relaxation 

of a qubit coupled to the oscillator ( Serban et al., 2010 1 . Spectral measurements have 

been reported both in the classical ( [Stambaugh and Chan, 20066 Almog et al., 2007) 



^T he results on the switching rate s for an additively modulated classical oscillator near the critical 
point | |Dykman and Krivoglaz, 1980[ l extend to the quantum regime if one replaces kgT — > hwo{n + 
1/2). 
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and quantum regimes ( [Wilson et ai, 2010[ ); see also the chapters by Chan and Stam- 
baugh and by Wilson et al. 

The spectra of interest are described by functions 

{{K{t)Lm) = ^ dt,{[K{t + u) - {K{t + umm - {Hmy, (7.29) 

recall that ujm is ujf for additive and lopjl for parametric modulation; is close to 
Wo- We will consider spectra near resonance, with |w| « wq; the operators K and L will 
be the ladder operators a or a^. In particular, the peak in the spectrum of spontaneous 
radiation emission by the oscillator is determined by Re ((a^, a))^ with u) ~ — wo, as in 
the absence of periodic modulation, cf. Mandel and Wolf (1995). A physical example 
is radiation from a nonlinear cavity, with the oscillator being the cavity mode modu- 



lated by an incident electromagnetic field ( Drummond and Walls, 1980} or excited by 
modulating the boundary of the cavity (Wilson ei aL, 2010). 



An additional weak resonant force A' exp(— iwt)-!- c.c. causes the oscillator to vi- 
brate at frequencies cj and 2a;M — i^. The vibrations at frequency cj are described 
by the scaled susceptibility x(<^)i which determines the corresponding displacement 
{bq) = (A'/2a;M)x(iw) exp(-iwt)-|- c.c. 



xM=*[((a,at))^-((at,a))l, 



(7.30) 



In the absence of modulation eqn (|7.30|) goes over into the standard expression for the 
oscillator susceptibility. Function xi'^) includes an extra factor 2a;M compared to the 
notation we used previously ( [Dykman et ai, 1994[ ). 

Spectra of a modulated oscillator have two major contributions. One comes from 
motion in the vicinity of stable vibrational states, where the oscillator spends much of 
the time, and the other comes from fluctuation-induced transitions between the states. 
We will consider them separately. 

7.4.1 Fluctuations about a stable vibrational state and local response 

The contribution to the oscillator spectra from motion near a stable vibrational state 
should be analyzed differently depending on the dimensionless oscillator relaxation 
rate n. For k ^ X (see below), one can linearize and then solve quantum equations 
of motion (|7.26p near the position {Qa,Pa) of a given stable state, similar to the 
classical case ( [Dykman and Krivoglaz, 1979 Dykman et ai, 1994) and to how it is 
done in the coherent state representation ( [Drummond and Walls, 1980[ ) . Equivalently, 
one can linearize in {Q — Q^,, P — Pa) the drift term in the master equation for the 
density matrix in the Wigner representation (|7.19p . The resulting contribution to the 



power spectrum is (Serban et ai, 2010 Dykman et ai, 2011) 



Re((a,at))[f) 



^2 (n + 1) {i^ - ci) + , 



nc2 



(1^2 



-(w-wm)- (7.31) 



The parameters in eqn (|7.3ip are expressed in terms of the scaled squared vibration 
amplitude in the considered stable state = Qa + -f? ' 
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(t'a)par = 



= «^ + irl - Art 



-1, (Ci)add = l-2?i, 

(Cl)par = - 2r^, (C2)pa 



(C2) 



add 



(7.32) 



where subscripts "add" and "par" refer to additive and parametric modulation, re- 
spectively; in eqn (|7.3ip we use the factor F/k as the frequency scale, this factor is 
independent of the decay rate F. In the limit of small damping, k — > 0, z^a goes over 
into the dimensionless frequency z/q of vibrations about the attractor, cqn (j7.14[) . 

The contribution of small-amplitude fluctuations to Re {{a\a)) _^ for uj close to 
Wo is also given by eqn (|7.31|) in which one should interchange n ^ n + 1. This 
relation, together with eqn ()7.30|) . give the imaginary part of the susceptibility x{'^)- In 
evaluating these local contributions we replace in the definition of the power spectrum 

CM Ht)) ^ (aW)^'^ = (2A)-i/2(Q^ + iPJexp(-zcjM0- 

The susceptibility for not too small k has the same form in the quantum and clas- 
sical case. It can be obtained by adding to the linearized in Q — Qa, P — Pa eqns (|7.26p 
a weak-driving term cx A', which oscillates in the rotating frame at frequency uj — ojm- 
One then finds the solution that oscillates at this frequency. Since the equations are 
linear, the noise term, which is the only term that has a different form in the classical 
and quantum case, drops out on averaging, giving for k 3> A 



Xa(w) 



V = -(w - cjm)- 



(7.33) 



For underdamped vibrations about the stable state in the rotating frame, where 
K <^ (but still K ^ A), the spectra (|7.3ip and Im Xa(w) have peaks at frequency 
detuning uj — wm = -^{V j k)vq. There are generally two peaks on the opposite sides of 
cjm and they generally have different amplitudes. Their shape is close to Lorentzian, 
with halfwidth F. As the ratio v^^jn decrease the peaks start overlapping. 



Additive driving, power spectrum 




Additive driving, absorption spectrum 

5r ' ' 




Parametric driving, power spectrum 




((Li-a>f)/6ai 



(26j-a)/,-)a)/,/F 



Fig. 7.8 The contributions to the scaled power spectra ^a{(jj) = Re {a /T) {{a\ a}}^^^^ and 
the susceptibility, Xa(^) = Im ('^/r)Xa('^) from fluctuations about stable vibrational states 
for K ^ A. The scaled decay rate is k = 0.3. For the additively driven oscillator, curves 1 
and 2 refer to the small- and large-amplitude vibrational states; the parameters are P — 0.15 
and n = 0.3. Curves 1 to 3 for the parametrically driven oscillator refer to /i = —0.6 and 
n = 0,0.3,0.5. 



The spectra are illustrated in Fig. 17.81 As seen from the figure and eqn (|7.31|) . the 
emission spectrum ((a^,a))^^ is symmetric for h = 0; this holds for both additive 
( [Drummond and Walls, 1980] ) and parametrically modulated oscillators. The onset of 
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emission for n = is related to quantum heating. In terms of quantum optics, if the 
modulating field is electromagnetic radiation and one considers emission of photons by 
the modulated oscillator, the emission can be thought of as resulting from multi-wave 
parametric process. Both the emission intensity and the positions of the spectral peaks 
depend on the modulation strength in a complicated way. 

An interesting feature of the susceptibility seen from Fig. 17.81 is that Im can 
become negative ( [Dykman and Krivoglaz, 1979[ ) . In the corresponding frequency range 
an additional weak field is amplified by the strong field, which can be also considered 
as a parametric multi-wave process. The amplification occurs in spite the fact that the 
absorption coefficient integrated over the whole spectrum is positive: from eqn (|7.33p 
/ dux'^{oj) ~ TT. This sum rule holds because, for any Fock state \N) (cf. Fig. 17. 3p . an 
induced dipolar transition up in energy has a larger amplitude [oc {N + 1)^^^] than 
down in energy (cx iV^/^). For a modulated nonlinear oscillator, the spectral regions of 
absorption and amplification are separated; weak field amplification generally occurs 
for both additive and parametric modulation. 

7.4.2 Fine structure of the spectra and the interference of transitions 

Interesting quantum effects emerge in the spectra for small k, where not only k <C I'o, 
but also At < A ^ I'o- In the whole range k <^ the stable states of the oscillator 
are at the extrema of g{Q,P). However, the approach of the previous section does 
not apply for k < A. To consider the range k ^ X, one has to take into account the 
nonequidistance of the quasienergy levels 5„, which results from the nonlinearity of 
vibrations about the extrema. For a given extremum at Q = Qo, -P = 0, to the leading 
order in the nonlinearity 



.90 + [XMn + 1/2) + X^Vn{n + l)/2] sgnggg. 



(7.34) 



where go — g{Qo,0). Parameter V — h'(){di> / dg) is determined by the slope of the 
dimensionless frequency v{g) of vibrations in the rotating frame with given g; the 
quantum correction ^ A to the frequency i>o — i'{go) is assumed to be incorporated. 
The parameters vq and V are plotted in Fig. 17.91 For small A the vibration nonlinearity 
is small, X\V\ ^ lyo. The quasienergy spectrum (|7.34|) is sketched in Fig. l7.10T a). 



Additive driving, small-amplitude stale 




Additive driving, large-amplitude state 



Parametric driving 




Fig. 7.9 The scaled eigenfrequency of vibrations about the stable states in the rotating frame 
I/O and the parameter of nonequidistance of quasienergy levels V, eqn (|7.34|) . for additive and 
parametric modulation. 

The power spectrum of the oscillator in a given stable state is formed by transi- 
tions |n) — >■ |n± 1) between neighboring quasienergy states in Fig. I7.10r a). One might 



Power spectra of modulated quantum oscillators 23 



expect that the spectrum is then a superposition of partial spectra that correspond to 
individual transitions, with width determined by the reciprocal lifetime of the respec- 
tive quasienergy states. However, to spectrally resolve the transitions one has to wait 
for time > (K/r)|Ay|~^. If this time becomes comparable or smaller than the lifetime, 
the transitions arc not independent, the transition amplitudes interfere. The overall 
oscillator spectrum is then formed by many interfering transitions. 

The typical number of states that contribute to the spectrum is determined by 
the effective Planck number fie, eqn (|7.15p . We assume that A|y|fie <^ z^o, so that 
all dimensionless transition frequencies are close to ii/g. Therefore the spectrum has 
well-separated peaks at w — wm ~ ±(r/K)i'o- Near the peaks, eqn (|7.29|) can be sim- 
plified using the interrelation (|7.14|) between operators Q,P and operators h,h^ that 
describe vibrations about the stable state in the rotating frame. The operators 6, W 
are the ladder operators of an auxiliary oscillator in thermal equilibrium, with scaled 
eigenfrequency i^O: energy spectrum (|7.34p . and temperature Te- The power spectrum 
of this oscillator is 

$^j,t (i^) = Re / dte'^^birpiO)) ^e^'^^/'^'Rc dte-*^^(&^(r)6(0)). (7.35) 



Jo Jo 

Equation (j7.35p is written for \i> — t'osgn^ggl <C I'o where $bf,t (i^) has a narrow peak. 
The occurrence of such a peak can be understood by noticing that, if one disregards 
nonlincarity and decay, 6(r) = exp(— iz/qt sgngQQ)6(0), and then $f,tit(z^) becomes a 
(5-function. The nonlincarity and decay lead to broadening of the (5-function. 

From cqns (|7.14|) and (|7.29|) . near the peak of the power spectrum of the original 
modulated oscillator on the high-frequency (low-frequency, for sgnggg < 0) side of 
the forced- vibration frequency, lu « wm + (r/K)t'o sgnggg, 

^Re ((a,at))(f) « ^e^'^«/^^Rc {{a\a))^^l « cosh^ ^.$,,t (i'), (7.36) 

where, as in eqns (|7.31|) and (|7.33p . = (K/r)(cj — wm). 

Similarly, near the peak in the power spectrum on the low-frequency (high-frequency, 
for sgnggg < 0) side of wm, where v = (K/r)(aj — wm) is close to — z^o sgn^QQ, 

^e^^°/-^^Re ((a,a^))i^) « ^Re ((at,a))L^i « sinhV*$h6t (-^^). (7.37) 

Equations (|7.30p . (|7.36l ). and ()7.37p also describe resonant peaks in the oscillator 
absorption spectrum Im XtSy)- The shapes of the peaks are determined by hmction 
$j{,t(z^). The peak near frequency wm + (r/K)i'o sgnggg corresponds to absorption of 
the additional field, Im Xa('^) > 0, whereas the one near wm — (r/K)t'o sgnggg, with 
Im Xaiyi) < 0, corresponds to amplification and has a smaller area. 

Equations ()7.36|) and ()7.37p reduce the problem of the spectra of modulated oscil- 
lator for weak damping, k <C i'o, to calculating the power spectrum of an auxiliary 
equilibri um oscillator $bbt(t^) ([Dykman et ai, 2011 1 . This problem was discussed pre- 



viously ( [Dykman and Krivoglaz, 1984 1 . The result for the spectrum near its maximum 
can be presented in the form 

'^bbf i'^) = {ne + l)Re Y.n=i (7-38) 
(/)(n, v) = 4n(A - l)"-i(A + l)-("+i) [n{2'^n - 1) - i{v - vq sg 
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Fig. 7.10 (a) A sketch of the quasienergy spectrum and dimensionless transition frequencies 
of weakly nonlinear vibrations about the stable state in the rotating frame for qqq > 0. 
(b) The power spectrum of the auxiliary oscillator near resonant frequency uo as given by 
eqn (|7.38|) for a comparatively large ratio of the level nonequidistance to the decay rate, 
i} — 10 and for qqq > 0. The curves 1 to 3 refer to fie = 0.1, 0.5 and 1.5. 

where A = [1 + *^(2ne + 1)] and H = [l + 2ii?(2ne + 1) - i?^] (Re H > 0). 
The parameter i? ~ {XV/2k) sgnggg is determined by the interrelation between the 
nonequidistance of the transition frequencies in Fig. IT.lOr a) and the decay broadening 
K of the quasienergy levels,. 

Equation ()7.38p represents the spectrum as a sum of effective partial spectra 
Re (j){n, v) that can be provisionally associated with transitions |n — 1) — \n) between 
the quasienergy levels in Fig. IT.lOf a). Functions 4'[n,v) depend on two parameters, 
•d and fie- The form of (j){n,v) is particularly simple for a comparatively large level 
nonequidistance or small damping, X\V\ 3> n, in which case 

A,( ^ n exphAt^o(n- l)/re] l ql _ \ 1^1 ^ , /n oa\ 

nn, v) w ,_ -r -, rr, \iy\ = A— > 1. (7.39) 

Here, i^((7„_i) = {vq + XVn)sgiigQQ is the frequency of the |n — 1) — > \n) transition in 
Fig. I7.10f a) and k„ = K[2n(2ne + 1) — 1] is the total halfwidth of levels \n — 1) and 
|n). Therefore Re 0(n, v) has a conventional form of a partial spectrum. 

For |t9| ^ 1 the overall spectrum $btt(i') has a fine structure. The intensities of 
the individual lines (|7.39p immediately give the effective quantum temperature Te. 
However, the fine structure is pronounced only in a limited range of the effective 
Planck numbers fie- This is seen from Eq. (|7.39p . For fie <C 1 only (f){l,i') has an 
appreciable intensity while Re Kcjy^njv) <^ 1 for n > 1. On the other hand, for large 
fie the linewidth k„ becomes large and spectral lines with different n overlap, starting 
with large n. Not only are they overlapping, but their shape is also changed compared 
to eqn (|7.39[) due to the interference of transitions. The evolution of the fine structure 
with varying as given by eqn (j7.38p is illustrated in Fig I7.10r b) . 

As ji?! decreases all partial spectra start to overlap, and for ji?! < 1 they can no 
longer be identified. As a result of the interference of transitions, in the limit — > 
we have (pin^v) cx and the spectrum has the form of a single Lorentzian peak 

of dimensionless halfwidth k, <&bbt (i^) = {fie + l)'^ ['^^ + (f^ ^ t'osgnggg)^] ^. This 
expression, with account taken of eqn (j7.36|) . agrees with eqn (|7.31|) in the range 
A 1 1^1 <C K <C vq where both apply. Generally, because of the nonlinearity of the auxil- 
iary oscillator, the shape of the spectrum depends on fie even where there is no fine 
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structure. The spectrum is non-Lorentzian and displays a characteristic asymmetry 
for large fie, where |i?|ne > 1- This asymmetry is described by cqn (j7.38p and provides 
a way of determining quantum temperature. 

Interference of transitions occurs also in various quantum oscillators and oscillator- 
type systems in the absence of modulation, from localized vibrations in solids to large 
spins in strong magnetic fields and Joscphson-j unction based systems, to mention but 
a few. In all these systems the spectra strongly depend on the interrelation between the 
level nonequidistance and the decay rate. Interference of transitions is important also 
for classical vibrational systems with fluctuating frequency, like nano- and microme- 
chanical resonators with a fluctuating number and/or positions of attached molecules 
dVig and Kim, 1999] [Yang et ai, 20TT ). The spectra of such systems can also be asym- 



metric and display a fine structure ( Dykman et al., 2010[ ). On the formal side, these 



recent results indicate that, for different systems and physical mechanisms, interfer- 
ence of transitions can be described by linear equations for coupled partial spectra. 
These equations are convenient for numerical solution. 

7.4.3 Supernarrow spectral peaks 

Along with fluctuations about stable vibrational states, quantum noise leads to occa- 
sional interstate switching discussed in Sec. 17.31 Important manifestations of switching 
are additional peaks in the oscillator power spectrum and susceptibility. The peaks are 
centered at the forced-vibration frequency wm and are supernarrow in the sense that 
their width is much smaller than the oscillator decay rate F. 

To describe the peak in the power spectrum we note that the populations wi and 
W2 = 1 — wi of the stable vibrational states 1 and 2 satisfy the balance equation 

dw^/dt^-W{wi-w^), W = W^l^^ +W^l^\ wi=W^l^^lW, (7.40) 

iij) . . . _ . 

where M4w is the rate of switching from zth to jth state and wi is the mean state 
population. Equation (|7.40|) has the same form as in the classical case, except that the 
switching rates are determined by quantum fluctuations. 

Fluctuations of the state populations lead to fluctuations of the expectation values 
of the operators a(t),a^(t) averaged over time ~ F~^, which switch between their 
stable-states values ai{t),a\{t), i = 1,2; in the lab frame in the zth stable state ai{t) = 
aiexp{~iuJMt) with = (2A)"i/2((5a> +iPaJ Pykman et al., 2011D . From eqn (fTjoj) . 
similar to the case of a classical oscillator ( [Dykman et ai, 1994 ), for the contribution 



from interstate switching to the power spectrum we obtain 

Re((a,at))[;") «Re((at,a))^^") « |ai - asj V«>2 + (c^ - c^m)']. (7.41) 

The typical width of the spectral peak (|7.4ip is given by the total switching rate W, 
it is exponentially smaller than the decay rate, cf. eqn ()7.16p . The intensity (area) of 
the peak is determined by the factor wi'W2- For a parametrically modulated oscillator 
for fi < fiB2 the populations of the states are equal by symmetry and this factor is 
equal 1/4. A peak in the power spectrum related to interstate transitions was seen in 



the radiation from a parametrically modulated microwave cavity (Wilson et ai, 2010). 

For additive driving, on the other hand, the populations wi and iB2 are exponen- 
tially different, and thus wi'W2 is exponentially small everywhere except for a narrow 
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(12) (21) 

parameter range where the switching rates Wsw and Wsw are ahuost equaL This 
range corresponds to a smeared kinetic "phase transition" ( [Bonifacio and Lugiato, 1978{ 
[Dykman and Krivoglaz, 1979| [Lugiato, 1984D , with the stable states playing the role 
of coexisting phases in a thermodynamic system. The onset of the supernarrow peak 
(|7.4ip is an indicator of the transition (Dykman et ai, 1994[). For a classical oscillator, 



such peak was observed by Stambaugh and Chan (200661 



The susceptibility of the oscillator also displays a supernarrow peak. The analysis of 
this peak for additive driving is similar to that for a classical oscillator ( [Dykman and Krivoglaz, 1979[ 
[Dykman et ai, 1994[ ). One notices that an extra force A' cxp{—iujt)+ c.c. with fre- 
quency very close to the strong- force frequency, |w — cji?| <^ T, can be described 
by making the parameter /3 in cqn (|7.7p slowly dependent on time, (3 — >■ = 
(3 + {2/3/A){A' exp[— z(a; — uJF)t] + c.c.} + . . . [we do not consider terms at the mirror 
frequency 2u!p — w in /3{t)]. One can then think of the switching rates becoming para- 
metrically dependent on time via /3{t). From eqn (|7.16|) . the major part of this time 
dependence comes from the modulation of the activation energies Rai and Ra2 for 
switching from states 1 and 2, respectively. The state populations also become time 
dependent, which gives the switching-induced contribution to the susceptibility 

X(^")(^) ~ wiW2—-^ a,-a2 (2;3/A)V2 _ g r^^) . (7.42) 

W — l(UJ — Up) LOp — U!q 

Function x^™H'^) displays resonant structure in a frequency range determined by 
the switching rate, jo; — wfI ^ W. The amplitude of x''™^ (i^) is proportional to a large 
factor i?Ai.2/A 3> 1. It also contains factor 'WiW2-, which steeply depends on the dis- 
tance to the kinetic phase transition. For a classical oscillator, the switching-induced 
peak of the response has been seen in the experiment ( [Chan and Stambaugh, 2006 



Almog et ai, 2007). The susceptibility of a parametrically modulated oscillator also 
displays a supernarrow peak for lu = wm with amplitude cx Ra/^ ^ 1. It was consid- 
ered previously for a classical oscillator ( [Ryvkine and Dykman, 2006[ ). 

The amplitude of the supernarrow peak of the susceptibility increases with de- 
creasing temperature, as the ratio Ra/^ increases. The very onset of this peak for 
low temperatures, as well as the supernarrow peaks in the power spectrum, is due to 
quantum fluctuation-induced interstate switching. 

7.5 Nonresonant modulation: oscillator heating and cooling 

We now briefly discuss oscillator dynamics in the presence of a moderately strong 
nonresonant modulation, where the modulation frequency is significantly different 
from the oscillator eigenfrequency wo- The dynamics of a modulated linear oscillator 
with linear in coordinate q and momentum p interaction with a thermal bath (|7.10|) was 
studied earlier ( [Schwinger, 1961 Zc ldovich et ai, 1970[ ). It can change dramatically if 



the interaction is nonlinear ( |Dykman, 1978 1 and/or the bath is also modulated 



A major effect of nonresonant modulation is the creation of a new relaxation mecha- 
nism. The relaxation discussed earlier in this chapter, see Fig. 17.31 was due to oscillator 
transitions between its neighboring energy levels. In a transition, the energy « Hujq is 
transferred to or taken from excitations in the bath. In the absence of modulation this 
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imposes on the oscillator the thermal distribution of bath excitations with energy fkoo, 
leading to the Boltzmann distribution pnn oc exjp{— N Hloq / ksT) over Fock states |A'^). 

In the presence of nonresonant modulation, in a bath-induced interlevcl transition 
the part huip of the energy can come from the modulation, see Fig. 17.11] The energy 
of the involved bath excitations is then h\ujQ ±ujf\- These excitations should impose 
on the oscillator their thermal distribution. It corresponds to oscillator temperature 
T* = Tojo/{oJo ± ujp). This temperature is positive if transitions \N) — >• — 1) 
correspond to emission of bath excitations and are thus more probable than transitions 
I A^ — 1) — >■ |iV), see Fig. l7.1lT a). (b). If, on the other hand, bath excitations are emitted 
in transitions \N — 1) — > |7V), then T* becomes negative. This case is sketched in 
Fig. l7.1lT c). It corresponds to relaxation processes in which the energy of a "photon" 
of the modulation goes into excitation of the oscillator and the bath. 



(a) 



(Op 
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Fig. 7.11 A sketch of modulation-induced relaxation processes leading to cooling (a), heating 
(b), and population inversion (c) in the oscillator; too, lof, and toh are the oscillator frequency, 
the modulation frequency, and the frequency of the bath excitation, respectively. A sufficiently 
strong modulation imposes on the oscillator the probability distribution of the bath modes 
involved in the scattering in (a) and (b) and leads to population inversion of low-lying energy 
levels and self-sustained vibrations of the oscillator in (c). 



Interestingly, the distribution of a nonresonantly modulated quantum oscillator 
over Fock states can remain of the Boltzmann form in a broad parameter range 
( [Dykman, 1978[ ). We will show this first for the case where the Hamiltonian of the 
interaction with the bath along with the linear in a, a' term (j7.10p has a quadratic 

Hl''>^q^h^^\ (7.43) 
The modulation term in the Hamiltonian of the isolated oscillator (|7.1I) is 

Hp = Hi^dd^ -qAcosupt; wf, |i.^_f - wo| > F, |7|(q^)/iuo- (7.44) 

We now go to the interaction representation using canonical transformation Unvit) = 
Tt exp{-ih-^ f dt'[Ho{t') + Hi,]}, where Ho{t) is given by eqns ([LT]) and ^A^, H^, 
is the Hamiltonian of the thermal bath, and Tt is the chronological ordering operator. 
For nonresonant modulation 

Ul{t)qUn,{t) w qo{t) + Aosc cosujFt, Aosc = A/{ujI - ujj,). (7.45) 



*The interaction 117.431 1 leads to nonlinear friction l |Dykman and Krivoglaz, 1984[ |, which can play 
an important role in oscillator dynamics, see the chapter by Moser et ah, and in particular determine 
the characteristics of self-sustained vibrations; we do not discuss this dissipation mechanism here. 
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Here, qoit) is operator q in the interaction representation in the absence of the modu- 
lating field. Similarly, operator a in the interaction representation is a sum of operator 
ao(i) calculated ioi A ~ and the terms oscillating as exp{±iujFt); the ac-Stark shift 
of the eigenfrequency ~ |7|^osc^/wo is assumed small compared to ujq. 

From eqns (|7.43p and (j7.45p . the nonlinear coupling to the bath in the interac- 
tion representation contains a term 2qo(t)h\^ (t)Aosc cosw^i, which has the same struc- 
ture as the linear coupling, eqn (|7.10p . except for the time-dependent factor cx Aqsc- 
Therefore the contribution of the field-induced scattering to the decay rate F is de- 
scribed by an expression similar to equation (|7.10p except that one should replace 
loq ujQ ± LOp- This gives the overall rate as V p = F + r+ + F_ — Fj, with 
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3 

dt{[h^\t)M^\mhe'^^''^^^^' K±WF>0); (7.46) 



rate — Fi is formally given by eqn (j7.46p for F_ in which — Wi? < 0; therefore Fi > 
and F_ = for ijjp — cjq > 0, whereas F_ > and Fi = for wq — lof > 0. The terms 
F± lead to the increase of the overall decay rate, whereas Fi reduces the rate and, if 
dominates, can make Tp negative. 

The mapping on the effectively linear in a, coupling suggests that, in the presence 
of nonrcsonant modulation, the master equation in the interaction representation has 
the conventional form for an oscillator with no modulation, cf. eqn (|7.1ip . 

dtp = -rp{np + l){alaop - 2aopal + pajao) - Tpfip{aoalp - 2alpao + paoaj,), 
rip = {Fn(cjo) + F+n(cjo + i^_f) + F_n(wo — ^f) + Ti [n{ujp - cjo) + 1]} /^F (7.47) 

where ap = ao{t), aj = aQ{t). Each relaxation mechanism in Fig. 17. Ill gives an additive 
contribution to eqn ()7.47|) . with the Planck number at the energy transferred to the 
bath in the corresponding elementary process. 

The stationary solution of eqn (|7.47p has the form of the Boltzmann distribution 
with temperature 

T* = hojo/kB \n[{np + l)/np]. 

This expression goes into the discussed above limits if one of the scattering mechanisms 
dominates, which requires that parameters F±,Fi significantly differ from each other. 

The difference between parameters F± , Fi is determined by the difference between 
the density of states of the thermal bath at the appropriate frequencies. It can be large 
if the oscillator is coupled to the thermal bath via an underdamped vibrational mode, 
with h'^^ being proportional to the coordinate of this mode qm ( [Dykman, 1978 ). The 



mode decay rate Tm should largely exceed F, |Ff|, but still it can be smaU compared to 
the mode frequency Um and jwm — wo|. One can then selectively tune cop to resonance 
and for example, for ujp fa LUm — 3> 1^0 achieve significant oscillator cooling. 

Similar effects occur and similar description applies if the coupling to the bath is 
linear in a, a' , eqn (|7.10l) . but nonrcsonant modulation is parametric, with Hamiltonian 
-f^par, eqn ()7.2|) . and ujp,\ujp — 2ujo\ S> F, |7|((7^)/wo- Parameters F±,Fi in this case are 
also quadratic in the modulation amplitude F and are determined by the correlators 
of /lb, /i], at frequencies ujQztujp. 
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The possibility of reducing the vibration temperature by modulation has been 
attracting much attention recently, in particular in the context of optomechanics, 
( Braginsky and Vyatchanin, 2002[ |Kippenberg and Vahala, 2008 ) . In an optomechan- 



ical system, the oscillator (a vibrating mirror) is coupled to a cavity mode, which 
is driven by external radiation. A quantum theory of cooling of the mirror in this 



case was developed by Wilson- Rae et al. (2007) and Marquardt et al. (2007). If the 



incident radiation is classical, in the appropriately scaled variables the coupling and 
modulation are described by Hamiltonians H'f^^ and respectively, with 

H'r^=c^qql^, H^;""^ ^ -q„AcosiJFt, (7.48) 

where q and qm are the coordinates of the mirror and the mode. In cavity optome- 
chanics one usually writes ifl™^ = Cm<l(4nam', the following discussion immediately 
extends to this form of the interaction. 

In the absence of coupling to the mirror the cavity mode is a linear system, hence 
<lm{t) — qmo(t) + [xm(wi?) exp(— iwFi) + c.c.JA/2, where qmo{t) is the mode coor- 
dinate in the absence of modulation and Xm{<^) is the susceptibility of the mode 
( Marquardt et a/., 2007[ ). The coupling in the interaction representation then 



has a cross-term oc qo{t)q^Q{t) exp(±iaJi?i). Since the cavity mode serves as a thermal 

(2) 

bath for the mirror, this term is fully analogous to the similar cross-term in Hi that 

f2l 

comes from modulation of the oscillator, with Crn(lmo{i) playing the role of hf^ (t). 
One can then describe the dynamics of the mirror by eqns ()7.46|) and ()7.47|) in which 

hl^\t) is replaced with Cmqmo{t) and A^^^ is replaced with \ Axm{uJF)\'^ ■ 

Another potentially important contribution to nonresonant cooling, heating, or 
excitation of the oscillator can come from the direct nonlinear interaction of the oscil- 
lator and the thermal bath with the modulation. For an electromagnetic modulation 
such interaction is due to nonlinear polarizability. For modulation Acosuipt the cor- 
responding interaction Hamiltonian is 

= -qh[f"> A cos u Ft. (7.49) 

For example, here can be the coordinate of a comparatively quickly decaying 
mode of a nanomechanical resonator. 

The effect of interaction ()7.49|) is again described by eqns (|7.46p and (|7.47|) in which 
h^\t) is replaced with h\f\t) and A^^^ is replaced with A^/A. If more than one of the 
above mechanisms is relevant, in calculating the rates of modulation-induced decay 
one should take into account the interference terms, i.e., write the overall effective 
coupling in the interaction representation as ao(t) exp{±iu;Ft)hh(t) -\- H.c. and then 
express the rates r± , Fi in terms of the commutator of hh (t) similar to eqn (|7.46|) . 

Nonresonant modulation modifies the power spectrum of the oscillator ( [Dykman, 1978[ 
|Dykman and Krivoglaz, 1984[ ) . Cooling of a nonlinear oscillator can lead to narrowing 
of its spectrum at frequency uiq even though F^? > F. Unless there are symmetry con- 
straints, the modulation can also lead to spontaneous emission of photons or phonons 
at frequencies \ujq ± ujf\- It can also lead to amplification of a weak external field at 
frequency ujf ~ for T* > 0. 
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7.6 Conclusions 

We described the dynamics of a modulated nonlinear oscillator for resonant and non- 
resonant modulation. Two types of resonant modulation were considered, an additive 
force with frequency ujp close to the oscillator eigenfrequency ujq and parametric mod- 
ulation with frequency lof close to 2wo- For resonant modulation, quantum dynamics 
in the rotating frame is characterized by three parameters, the scaled modulation in- 
tensity, decay rate, and Planck constant, which are given in Table 17.11 Of primary 
interest was the parameter range where the oscillator displays bistability of forced 
vibrations. Fluctuations were assumed small on average, so that the smearing of the 
classically stable states in phase space is small compared to the interstate distance. 
For T > the system lacks detailed balance 

Relaxation of the modulated quantum oscillator is accompanied by quantum noise. 
It leads to a finite- width distribution over quasienergy states even for T — >• 0. We found 
the distribution for weak damping, where the width of the quasienergy levels is small 
compared to the level spacing. Near its maximum the distribution is Boltzmann-like. 
The far tail is of non-Boltzmann form. It determines the exponent in the rate of switch- 
ing between the stable vibrational states. The switching occurs via transitions over the 
effective barrier that separates the states in phase space. We called this quantum ac- 
tivation and studied the corresponding effective activation energy. Remarkably, even 
for T — )■ it is smaller than the exponent for switching via tunneling. Therefore inter- 
state switching occurs via quantum activation, not tunneling, unless the decay rate is 
exponentially small. We also found, using a different method, the switching rate close 
to bifurcation points. Its scaling with the parameters is given in Table [7^ 

We found the power spectra and the susceptibility of resonantly modulated oscil- 
lators. The spectra have a characteristic shape, which depends on the interrelation 
between the decay rate, the quasienergy level spacing, and the nonequidistance of the 
quasienergy levels. Where the nonequidistance exceeds the level width, the spectra 
display a fine structure, which sensitively depends on the effective temperature of the 
quasienergy distribution near its maximum. They can also display a characteristic su- 
pernarrow peak where the stationary populations of the coexisting vibrational states 
are close to each other. 

An interesting effect of nonresonant modulation of the oscillator is that it can sig- 
nificantly change the oscillator distribution over the Fock states, leading to heating, 
cooling, or excitation of self-sustained vibrations depending on the modulation fre- 
quency and the coupling to the thermal bath. We show that different coupling and 
modulation mechanisms can be described in a similar way and that the distribution 
can remain of the Boltzmann form in a broad range of oscillator energies. 

This research was supported in part by DARPA through the DEFYS program and 
by the NSF, Grant EMT/QIS 089854. 
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